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random impurity potential, respectively. Within both models, we find that the anomalous Hall con- 
ductivity has a finite real value at zero frequency, exhibits singularities at the threshold of photon 
absorption across the superconducting gap, and decays as some power of the high frequency uj. The 
Hall conductivity increases linearly with the decrease of temperature below the superconducting 
transition and saturates at zero temperature. Using our results for the high-frequency Hall conduc- 
tivity, we estimate the polar Kerr angle for light reflection from the material and compare it with 
the experimental measurements in Sr2Ru04 by Xia et al., Phys. Rev. Lett. 97, f67002 (2006). 
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I. INTRODUCTION 

Recent development of an optical apparatus based on 
the Sagnac interferometer by the group of Kapitulnik 
led to spectacular discoveries of the spontaneous time- 
reversal symmetry breaking (TRSB) in the unconven- 
tional superconductors Sr 2 Ru04 [1] and YBa 2 Cu306+ x 
[2] . These experiments measured the rotation angle 9 k of 
the polarization plane of light upon normal reflected from 
the material surface. A non-zero value of the Kerr angle 
9k indicates the TRSB in the sample [3]. Positive and 
negative values of 9k represent clockwise and counter- 
clockwise rotation of polarization and indicate two pos- 
sible ways of spontaneous breaking of the time-reversal 
symmetry. In S^RuCi, a non-zero 9k appears at the 
superconducting transition temperature T c — 1.5K, so 
the TRSB in this material is clearly related to the onset 
of superconductivity [1]. In contrast, in the underdoped 
YBa2Cu306+ x , a non-zero 9k appears in the tempera- 
ture range where the pseudogap develops, which is much 
higher than the superconducting transition temperature 
[2]. Thus, the TRSB in cuprates appears to be due to 
a second-order phase transition with an order parameter 
unrelated to superconductivity. 

In this paper, we focus on manifestations of the TRSB 
in Sr2RuC>4, which is a quasi-two-dimensional (Q2D) 
metal consisting of weakly coupled conducting layers in 
the (x,y) plane [4, 5]. It was proposed theoretically 
that the superconducting pairing in Sr2Ru04 is spin- 
triplet with the chiral p x ± ip y orbital symmetry [6] . In 
this state, the Cooper pairs have the angular momentum 
L z = ±1 pointing perpendicular to the conducting lay- 
ers of Si'2Ru04. This order parameter breaks the time- 
reversal symmetry and is the two-dimensional analogue 
of the A phase in superfluid 3 He [7-9]. There is sub- 
stantial experimental evidence indicating that Sr2Ru04 
is indeed a spin-triplet [10, 11] p-wave [12] superconduc- 



tor. However, the chiral orbital symmetry of the order 
parameter and the violation of the time-reversal sym- 
metry are still under debate [13]. The early evidence 
for the TRSB came from the muon spin-relaxation mea- 
surements [14]. Recent observations of the spontaneous 
polar Kerr effect in optics [1] and dynamical domains 
of the superconducting order parameter in Josephson 
junctions [15] give strong evidence for the chiral pairing 
with the TRSB [16, 17]. However, the scanning SQUID 
and Hall probe experiments [18, 19], designed to search 
for domains with opposite chiralities on the surface of 
Sr2RuC>4, did not find any evidence for such domains 
and for the TRSB. The discrepancy between these ex- 
perimental results shows that the TRSB in superconduc- 
tors is not fully understood and requires further investi- 
gation [13]. Chiral superconductors with the TRSB are 
expected to have many unusual properties, e.g., certain 
vortex excitations with zero-energy Majorana modes in 
the core [20, 21]. These non-local quasiparticle excita- 
tions obey non-Abelian statistics and are studied in the 
context of topological quantum computation [22]. The 
fractional quantum Hall state at the filling factor 5/2 
was proposed to be analogous to the p x + ip v supercon- 
ductivity [20, 23]. 

Here we study an important manifestation of the 
TRSB in the superconducting state: the emergence of 
anomalous (spontaneous) Hall conductivity a xy in the 
absence of an external magnetic field. This study is mo- 
tivated by the experimental measurement [1] of the polar 
Kerr effect in Sr 2 Ru04. According to the textbook calcu- 
lations using classical electrodynamics [24], the Kerr an- 
gle 9k is determined by the ac Hall conductivity cr xy (uj), 
where u) is the frequency of light in the experiment. Thus, 
the experimental observation [1] of a non-zero 9k at 
T < T c signifies appearance of a non-zero anomalous ac 
Hall conductivity in the material. However, a theoretical 
calculation of the anomalous ac Hall conductivity for a 
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chiral p x + ip y superconductor turned out to be a rather 
non-trivial problem and led to some controversies briefly 
reviewed in the next Subsection. 



A. Theories of the anomalous ac Hall conductivity 
in a clean chiral p x + ip y superconductor 

Most of the previous theoretical calculations of the 
anomalous Hall conductivity focused on the clean limit 
for a chiral translationally-invariant superconductor in 
the absence of impurities. However, a straightforward 
calculation of the current-current correlation function 
shown in Fig. 1 and involving interaction with the A x 
and A v components of the electromagnetic vector poten- 
tial gives zero result [7, 25-32]. This implies that the Hall 
conductivity <r xy in a clean chiral superconductor is zero, 
even though a non-zero value is nominally permitted by 
the TRSB. Kim et al. [33] obtained a non-zero value of 
<j xy for an anisotropic superconducting pairing, but their 
calculated conductivity tensor is symmetric a xy — a yx . 
So, it does not represent the antisymmetric Hall conduc- 
tivity (a xy = —cr yx ) and does not contribute to the polar 
Kerr effect. 

On the other hand, it was found that the charge- 
current correlation function, represented by the Feynman 
diagram with the scalar electromagnetic potential A and 
the vector potential A x or A y , is non-zero in chiral super- 
conductors [7, 26-32]. It describes the magneto-electric 
effect, i.e., a change in the electric charge density in re- 
sponse to an applied magnetic field B z . As a consequence 
of this effect, a magnetic vortex would acquire the frac- 
tional electric charge e/4 [34]. A setup for an experi- 
mental detection of the magneto-electric effect in chiral 
superconductors was proposed in Ref . [32] . 

In conventional (non-superconducting) quantum Hall 
systems, the magneto-electric effect is directly related to 
the Hall conductivity via the Streda formula [35]. How- 
ever, in superconductors, this relation is invalid, and the 
magneto-electric effect does not imply the Hall effect. In- 
correct claims were made in Rcfs. [7, 31, 36] on the basis 
of the magneto-electric effect that clean chiral supercon- 
ductors have a non-zero anomalous Hall effect. It was 
shown in Refs. [26-29, 32, 37] that, when the contribu- 
tion from the collective current of the superconducting 
condensate is properly taken into account in a gauge- 
invariant manner, the Hall conductivity vanishes. This 
result is also in qualitative agreement with the conclu- 
sions of Ref. [38]. In principle, a non-zero value of the 
Hall conductivity a xy (q 7 uj) can be obtained for a non- 
zero in-plane wave vector q [26-29, 32]. However, the 
estimates done in Ref. [32] using q derived from the laser 
beam diameter in the experiment [1] give the value of 
<j xy many orders of magnitude smaller than what is nec- 
essary to explain the experiment. A non-zero, but very 
small Kerr angle was obtained in Ref. [39] by considering 
the orbital collective modes in a p x + ip y superconduc- 
tor. However, the estimated value of Ok is too small to 



explain the experiment [1]. 

Refs. [40, 41] claimed that, in unconventional super- 
conductors, it is necessary to make the Peierls-Onsager- 
style substitution of the momentum p — > p — eA in the 
momentum-dependent superconducting order parameter 
A(p). As a result, a new vertex of interaction with the 
electromagnetic field A would appear in the Hamiltonian 
of the superconducting system, and this vertex may give 
a non-zero value to a xy . However, this substitution is 
wrong for a number of reasons, as explained in detail 
in Appendix A. The absence of such substitution was 
shown in several papers, including Refs. [7, 23, 42, 43] 
and was recently recognized in Ref. [44]. 

In fact, the vanishing of the anomalous Hall conductiv- 
ity for a clean superconductor can be understood on gen- 
eral grounds using Galilean invariance of the system [20] . 
In the absence of an external magnetic field, the electric 
field E x applied in the x direction cannot induce a center- 
of-mass motion of the electron gas in the transverse y 
direction no matter whether the pairing between elec- 
trons is chiral or not, because there is no total external 
Lorentz force acting in the y direction. Thus, the experi- 
mental observation of the polar Kerr effect [1] cannot be 
explained within any model of a clean Galilean-invariant 
chiral superconductor [45]. 

The general argument given above uses the Galilean 
invariance of a translationally-invariant system. How- 
ever, electrons in a crystal are subject to a periodic lat- 
tice potential, which breaks translational symmetry and, 
thus, may invalidate the general argument. In fact, there 
are well-known examples of periodic systems, such as the 
TRSB topological insulators and metals, where a topo- 
logically non-trivial band structure produces the anoma- 
lous Hall effect [46, 47]. However, this effect has nothing 
to do with superconductivity and, thus, is not relevant for 
the Kerr effect in Sr 2 RuC>4. Relation between chiral su- 
perconductors and the TRSB topological insulators and 
metals is discussed in more detail in Appendix B. 



B. The role of impurities in the anomalous Hall 
effect 

Following the general argument presented in the pre- 
ceding Subsection, we realize that, in order to obtain a 
non-zero anomalous Hall conductivity, it is necessary to 
identify a physical mechanism for producing an external 
force on the electron gas in the transverse y direction. 
Scattering on impurities breaks translational invariance 
and provides a mechanism for momentum transfer be- 
tween electrons and the crystal lattice, which may gen- 
erate a net force in the y direction, if the pairing be- 
tween electrons is chiral. Making an analogy between a 
Cooper pair with the angular momentum L z = 1 and 
a spinning baseball, we observe that the baseball does 
not deflect sideways when flying in vacuum, but does de- 
flect in the air, because friction with the air generates the 
transverse Magnus force. Similarly, scattering on impuri- 
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ties can provide an effective friction between the spinning 
Cooper pairs and the crystal lattice and, thus, generate a 
transverse force. It should be emphasized that a periodic 
lattice potential alone, in the absence of impurities, docs 
not produce a non-zero Hall conductivity for a chiral p- 
wave superconductor [32] . It is the random scattering on 
impurities which provides a mechanism for momentum 
relaxation in the electron gas. 

The effect of disorder on the anomalous Hall con- 
ductivity in a p x + ip y superconductor was studied by 
Goryo [48]. He considered the skew-scattering processes, 
which involve triple scattering events on a given impu- 
rity. This process requires existence of a non-zero third 
moment of the impurity potential averaged over random 
realizations. Since the third moment is zero for the com- 
monly used Gaussian distribution, we call this type of 
distribution the non-Gaussian disorder. A similar model 
of disorder was used in the context of the anomalous 
Hall effect [49-52], as well as the spin Hall effect [53-55], 
for non-supcrconducting materials. In fact, Sinitsyn [47] 
was the first to point out that impurity scattering may 
contribute to the anomalous Hall effect in chiral super- 
conductors. 

Goryo [48] calculated the anomalous ac Hall conduc- 
tivity <r xy (cj) for a p x + ip y superconductor in the limit 
of high frequencies u) 3> Ao, where 2Ao is the supercon- 
ducting gap, and for T close to T c . This limit is rele- 
vant for the experiment [1], which utilized u = 0.8 eV, 
whereas the gap at T — can be estimated from the BCS 
relation as 2Ao = 3.5 T c — 0.46 mcV. Goryo found that 
the anomalous Hall conductivity is predominantly real in 
this limit and behaves as e' xy (u)) oc 1/w 3 . However, this 
result is inconsistent with the general causality relation 
<t(u>) = <T*(-u>) [56]. 

In this paper, we calculate the real and imaginary 
parts of the antisymmetric conductivity tensor a xy (u) = 
—<r yx (uj) originating from skew-scattering on impurities 
for the full range of frequency u and temperature T con- 
sistently with the causality relation. We find that the 
factor of i = was overlooked in the calculation of 

Ref . [48] . Then we consider a more conventional model of 
Gaussian disorder and calculate the Feynman diagrams 
where electrons scatter on two different impurities. We 
show that these diagrams also give a non-zero anoma- 
lous Hall conductivity for a chiral superconductor if the 
particle- hole asymmetry is taken into account. The role 
of the particle-hole asymmetry was discussed in Ref. [39] 
for chiral p-wave superconductors and in Refs. [57-60] 
for the Hall effect anomaly associated with vortex mo- 
tion in the high-T c superconductors. Then we compare 
the expressions for a xy (u>) obtained for these two models 
of Gaussian and non-Gaussian disorder and discuss the 
dominant contribution at the high frequencies relevant 
for the Kerr effect measurements [1]. 

The paper is organized as follows. In Sec. II, we first 
show that a xy vanishes for a clean chiral p x + ip y super- 
conductor. Then we show that the lowest (second) order 
diagrams in the strength of the impurity scattering also 



give a vanishing a xy . After that, we present a general 
discussion of the higher-order contributions, which lead 
to a non-zero a xy . In Sec. Ill, we calculate a x 3 y \u>) orig- 
inating from the skew-scattering diagrams for the non- 
Gaussian model of disorder. In Sec. IV, we calculate 
a x 4 y (<j) for the Gaussian model of disorder in the pres- 
ence of particle-hole asymmetry. The implications of our 
results for the experiments in Sr2Ru04 are discussed in 
Sec. V, and conclusions are given in Sec. VI. Technical 
details of calculations are relegated to the Appendixes. 



II. THEORETICAL MODEL AND GENERAL 
DISCUSSION 

In this Section, we first introduce a theoretical model 
for calculation of the anomalous ac Hall conductivity. 
Then, after a brief discussion of the clean limit, we 
present a general discussion of the effects of impurity 
scattering on the anomalous Hall conductivity. 

We use the electromagnetic gauge where the scalar po- 
tential A is set to zero, and the electromagnetic field is 
characterized by the transverse component of the vector 
potential A. Within the linear response approach, the 
current j appearing in response to an infinitesimal vec- 
tor potential A is 

j(q,u) =Q(q,w)A(q,w), (1) 

where Q{q 1 uo) is the electromagnetic response tensor. 
From this tensor, we can obtain the conductivity tensor 
aij (q, u>) relating j with the electric field E = —d t A 



In the experiment [1], the light beam is incident along the 
z axis perpendicular to the two-dimensional (2D) con- 
ducting layers of Sr 2 Ru04, and the vector potential A is 
polarized in the (x, y) plane of the layers. So, we calculate 
the off-diagonal response function Q xy (uj) = —Q yx (u>) at 
q = 0, which determines the ac Hall conductivity a xy (oj) 
and the Kerr angle Ok [24]. 



A. The model of a chiral p x + ip y superconductor 

The triplet superconducting pairing is characterized 
by the vector d, which determines spin polarization of 
the triplet Cooper pairs [61]. For Sr 2 Ru04, we consider 
the case where the vector d has a uniform, momentum- 
independent orientation [9]. Selecting the spin quanti- 
zation axis along the vector d, we obtain the represen- 
tation [32] where the triplet Cooper pairing takes place 
between electrons with the opposite spins [62]. For the 
orbital symmetry, we consider the chiral pairing potential 
A(p) = A (p x ± ip y )/pF, where p x , p y , and p F arc the 
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in-plane electron momenta and the Fermi momentum. It 
is convenient to write this pairing potential in the form 



and set A x 
calculations. 



A(p) 



A X P X + lAyPy 



(3) 



± y — A /pf only at the end of the 
The sign of the product 



sign(A a; A !y ) 



(4) 



and A y is different 



represents the sign of the order-parameter chirality. No 
tice that the dimensionality of A 
from that of A . 

In the Matsubara representation, the action 5 e i of 
the electron system can be written [32] as a 2 x 
2 Nambu matrix acting on the spinor i/>(p, w{) — 
[ip<t(p,vvi),ip^(—p,zui)], where ip and ipi are the destruc- 
tion and creation operators of the electrons with the mo- 
mentum p, the fermionic Matsubara frequency w\, and 
the spin projection | or J. 



S c i = ^2 V^P, ™i) [i^iTo - £(p) r 3 



P X A X fl +PyAyT 2 } tp(p, Wl). 



(5) 



Here 71,2,3,0 are the Pauli matrices and the unity matrix 
acting on the spinor ip. We set the Planck and Boltz- 
mann constants to unity: T% = 1 and ks = 1- The 
function £(p) represents the electron energy dispersion 
counted from the chemical potential /i. Interaction of 
electrons with the electromagnetic field is described by 
the action S e m 



S cm = ~e ^2 v{p)-A(iu n )i)} ] (p, w\ +w n )f tp{p,wi) 



(6) 



Here e is the electron charge, v(p) — d£ ) (p)/dp is the 
electron velocity uj n is the bosonic Matsubara frequency 
of the electromagnetic field, and the wave vector q of the 
electromagnetic field is set to zero. 

Sr 2 RuC>4 is a Q2D metal consisting of weakly coupled 
conducting layers. Its electron dispersion £(p) can be 
written as 



£(p) = £\\(Px,P y ) - 2t ± cos(p z d) - /i. 



(7) 



The first term in Eq. (7) represents the in-plane electron 
dispersion, and the second term the tight-binding out- 
of-plane dispersion with the tunneling amplitude t± and 
the interlayer distance d. Generally, the Fermi surface 
of Sr 2 Ru0 4 is rather complicated and consists of three 
sheets [5] . To simplify presentation, we consider only one 
sheet and assume that the in-plane dispersion is isotropic, 

i.e., £|| depends only on p = ^Jp x + p yl which is a good 

approximation for the 7 sheet [5] . 

We assume that the interlayer tunneling is weak tj_ <C 
/i, so the Fermi surface is a slightly corrugated cylinder 
extended along the p z direction with the average radius 




FIG. 1: The Feynman diagram representing the anomalous 
Hall conductivity of a clean chiral p-wave superconductor. 
The wavy lines represent an external electromagnetic field, 
and the solid lines are the electron Green's functions in the 
Nambu representation. 



Pf in the plane [9]. For a given value of the out-of-plane 
momentum p z , Eq. (7) shows that the system can be 
treated as a 2D metal with the effective chemical po- 
tential fb(p z ) = M + 2t±cos(p z d) weakly dependent on 
p z . Such an effective 2D description is possible, because 
the pairing potential A{p) in Eq. (3) does not depend 
on p z , and the vector potential A in Eq. (6) is polar- 
ized in the (x, y) plane. So, in the rest of the paper, we 
calculate the anomalous Hall conductivity a xy per one 
layer for a purely 2D electron system with the disper- 
sion £(p) = s\\(p) — p. A generalization to the Q2D case 
involves trivial additional averaging over p z within each 
fermion loop of the Feynman diagrams. The bulk Hall 
conductivity is obtained by dividing the 2D result for a xy 
by the interlayer distance d. 



B. Anomalous Hall conductivity of a chiral 
superconductor in the clean limit 

The electromagnetic response function Q xy (I) for the 
model (5) and (6) is given by the Feynman diagram in 
Fig. 1. The analytic expression for this diagram is 



Qxyi^n) = e 2 Tr v x {p) G{zui,p) v y (p) G(zui + uj n ,p) 



(8) 



where Tr denotes the trace over the Nambu space, as well 
as the sum over the internal momenta p and fermionic fre- 
quencies wi. The Green's function for the chiral p-wave 
superconductor in the Nambu representation is obtained 
by inverting the kernel in Eq. (5) 



G{w h p) 



iWlTQ + $pf 3 + P X A X T\ - PyA y T 2 

w?+e P + Px A x + Py ]Ai ■ 



(9) 



Here we write the argument p of the dispersion £(p) as 
the subscript in order to shorten the mathematical equa- 
tion. 

Using Eq. (9) for calculation of the trace over the Pauli 
matrices f in Eq. (8), we find 



f d? 

Q xy = 2e 2 T J T^VxVy 



iwi(iwi-\-iui r , 



e 2 



(10) 
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where E p = + Aq is the quasiparticle energy. It is 

clear that the response function Q xy in Eq. (10) vanishes 
upon integration over the orientation of p, because v x 
is an odd function of p x , and v y is an odd function of 
p y . Even if an anisotropic model with A x ^ A y is con- 
sidered [33], still a similar calculation gives a symmetric 
tensor Q xy — Q yxi which does not represent the Hall con- 
ductivity [25] and is not relevant for the experiment [1]. 

C. Disorder-induced anomalous Hall conductivity 

As discussed in Sec. IB, in order to obtain a non-zero 
anomalous Hall conductivity for a chiral superconductor, 
it is necessary to include the effect of disorder. Thus, we 
add the impurity scattering term to the action of the 
system 



*5*imp — 



E 

q,p^i 



V imp (q)ip r (p + q,wi)T 3 tp(p,zui), (11) 



where Vi mp (q) is the impurity potential written in the 
momentum representation. We assume that the domi- 
nant scattering mechanism comes from the short-range 
disorder. The first and the second moments of the prob- 
ability distribution function of the impurity potential 
Vimp(q) are 



(V imp (q)) = 0, 
(Vimp(q) V imp (q')} = rnulS(q + q'). 



(12) 



Here the averaging is performed over different realiza- 
tions of the disordered potential, and m and uq are the 
2D concentration of impurities and the strength of the 
disorder potential, respectively. 

It is well known that impurity scattering suppresses un- 
conventional superconducting pairing when the parame- 
ter A r (where r is the quasiparticle scattering time) 
becomes of the order of unity, because the Anderson the- 
orem does not apply to non-s-wave pairing [63]. Sup- 
pression of superconductivity by disorder has been ex- 
perimentally observed in Sr 2 Ru04 [64] and is one of the 
arguments in favor of the unconventional pairing sym- 
metry. Therefore, we study the case where the concen- 
tration of impurities is very low, so that A r » 1, and 
the superconducting pairing is not significantly affected 
by impurities. In fact, Sr2Ru04 is a very clean stoichio- 
metric material, where impurities are not introduced in- 
tentionally, and physical origin of residual disorder is not 
clear. Nevertheless, the presence of a small, but non-zero 
concentration of impurities can still induce a non-zero 
anomalous Hall effect. For such a low concentration of 
impurities, the anomalous Hall response can be studied 
perturbatively in the strength of the disorder potential 
and would be dominated by the the lowest-order non- 
vanishing diagrams. 

The lowest-order contributions to the anomalous Hall 
conductivity appear in the second order in Vi mp (g) as 




FIG. 2: The lowest-order diagrams in impurity scattering for 
the anomalous Hall conductivity. The dashed lines represent 
scattering potential from the same impurity, shown by the 



shown in Fig. 2. The diagram a) represents a self-energy 
correction to the Green's function of the electrons due 
to disorder. In p-wave superconductors, the momentum 
average of the gap is zero: J2 P A(p) = 0. Thus, the 
self-energy corrections modify only the To and T3 com- 
ponents of the Green's function (9). However, such a 
modification does not change the Nambu structure of the 
current-current correlation function (8). Therefore, sim- 
ilarly to the clean case, the diagram a) gives a vanishing 
contribution to the anomalous Hall conductivity. 

On the other hand, the diagram b) in Fig. 2 has a non- 
trivial structure and deserves a more detailed discussion. 
The analytical expression for the diagram b) is 



] x 2 y\uj n ) = mul Tr k x %K y f?, 



(13) 



where Tr is taken over the Pauli matrices and the internal 
frequency mi. The effective vertices A^ and A y are 

A x = -e^2,G(mi + w„,pi) v x (Pi) G(m h pi) (14) 
pi 

ep 2 F A r° d£ i(2zvi+u n ) 
= - — A xTl J 



iep 2 F A x , 
— a T i 



^ 2n \{wi + Lu n ) 2 + E 2 ][zuf + Ep] 
1 1 \ 



y/rf + A o V(zui+u; n ) 2 + A 2 



and 



= -e^2G(mi,p 2 ) v y (p 2 ) G(mi +ui n ,p 2 ) (15) 

P2 

.ep%. . f D d£ {2w l +uo n ) 
z J -1 



iep 2 F A y 



_ jJD 2n [( m +u; n ) 2 + E 2 ][tf + E 2 ] 
1 1 \ 



4^„ T2 ^ y/^f + A 2 + LU n ) 2 + A 2 ) ' 

In going from the first to the second lines in Eqs. (14) 
and (15), we integrated over the angular orientation of 
the electron momentum p. For the integration in the 
radial direction, we make the linearized approximation 
for the 2D electron dispersion 



£(p) w v f (p-Pf), 



(16) 



where vp is the Fermi velocity. This approximation is 
justified, because the integrals converge in the vicinity of 
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FIG. 3: The skew-scattering diagrams contributing to the 
anomalous Hall conductivity. The three connected dashed 
lines represent the third moment (17) of the impurity poten- 
tial. 

the Fermi surface. Then, the integration over dp can be 
replaced by integration over = vpdp. The limits of in- 
tegration over £ are given by the BCS cutoff energy ud, 
which is determined, presumably, by the energy scale of 
ferromagnetic fluctuations inducing the p-wave supercon- 
ductivity in Sr2Ru04. Throughout the paper, we assume 
that ujd 3> ui, A and neglect corrections to Eqs. (14) and 
(15) of the order 0(cu^ 3 ). 

Substituting Eqs. (14) and (15) into Eq. (13), we find 

that Q x 2 y^ — because Tr{fif3f2f 3 } = 0. So, there is 
no contribution to the Hall conductivity to the second 
order in V^ mp . Thus, we have to consider the higher-order 
diagrams in order to obtain a non-zero Hall conductivity. 

D. Anomalous Hall conductivity due to the 
higher-order diagrams in disorder potential 

In the rest of the paper, we study the higher-order con- 
tributions in disorder potential. This requires to specify 
the assumptions about the probability distribution of the 
disorder potential V] mp . Below, we consider two models 
with the Gaussian and non-Gaussian distributions of the 
disorder potential. In the case of the Gaussian distribu- 
tion of disorder, it is enough to specify the second-order 
cumulant given in Eq. (12). All other higher-order mo- 
ments are either zero (odd in V[ mp ) or can be expressed 
in terms of the second-order cumulant (even in V5 mp ). In 
the case of a non-Gaussian distribution of disorder, the 
odd cumulants may be non-zero. In particular, the third 
moment (skewness) within a non-Gaussian model is 

(^ m p(q) Mmp(t?') V imp {q")) = K 3 n lU 3 5{q + q' + q"). 

(17) 

Here K3 is a dimensionless parameter characterizing 
the skewness, which varies from to 1 depending on 
the deviation of the actual distribution function from 
the Gaussian. For example, K3 = 1 corresponds to 
randomly distributed impurities generating short-range, 
delta-function potentials of the equal strength u . The 
contributions to the Hall conductivity to the first order 
in the third moment (17) are shown in Fig. 3. 

The skew-scattering diagrams, similar to those shown 
in Fig. 3, were studied for the anomalous Hall effect in 
ferromagnets [50-52] and for the extrinsic spin Hall ef- 
fect [53-55] in semiconductors. In the case of a chiral 




zui +ui n m + u„ 



FIG. 4: The lowest-order diagrams contributing to the Hall 
conductivity within the Gaussian model of disorder with the 
particle-hole asymmetry. 

p-wave superconductor, the anomalous Hall conductivity 
due to the skew-scattering diagrams of Fig. 3 was first cal- 
culated by Goryo [48] for high frequencies u> 3> Ao and T 
close to T c . In Sec. Ill, we calculate the real and imagi- 
nary parts of a xy {uo) originating from the skew-scattering 
diagrams shown in Fig. 3 for the full range of frequencies 
lo and for any temperature T. 

For the Gaussian model of disorder, the higher-order 
diagrams contributing to the Hall conductivity are shown 
in Fig. 4. We calculate the real and imaginary parts of 
<7 xy (oj) originating from these diagrams in Sec. IV. We 
show that these diagrams give a non-zero contribution 
to the Hall conductivity only when we take into account 
the particle-hole asymmetry, which is discussed in more 
detail in the next Subsection. 



E. Sign of the anomalous Hall conductivity 

In this Subsection, we use the symmetries of the prob- 
lem to discuss what factors determine the sign of the 
anomalous Hall conductivity. 

First, let us discuss the time-reversal operation. The 
Hall conductivity changes sign upon the time reversal: 
a xy — ► — <j xy . On the other hand, the time-reversal oper- 
ation also results in complex conjugation of the Hamilto- 
nian. For a chiral p-wave superconductor, it means that 
A(p) — > A*(p), so the pairing potential (3) changes its 
chirality. Thus, we conclude that 

a xy (A) = -a xy (A*), (18) 

i.e., changing the sign of chirality of the superconducting 
order parameter changes the sign of the anomalous Hall 
conductivity, so <r xy oc s xy , see Eq. (4). 

When we employ the commonly used linearized ap- 
proximation (16) for the electron dispersion, the prob- 
lem acquires an additional symmetry upon exchange of 
the electron and hole operators, tp(p) — > tp^{p) and 
ip'(p) — > ip(p)- This operation preserves the fermion an- 
ticommutation relations, so it is a canonical transforma- 
tion. Upon the particle-hole transformation and subse- 
quent commutation of the tp and ift operators, the Hamil- 
tonian of the kinetic energy of electrons ^2 £p 1 /^(^ 5 )' ^ / ; (^ > ) 
transforms to the original form (up to an additive con- 
stant), but with the sign change of the energy disper- 
sion £ p — ► — £ p . If, in addition, we use the approxi- 
mation (16) and change the radial momentum variables 



7 



p — > 2pi? — p in the sum, the kinetic energy returns to the 
original form, i.e., becomes invariant under the particle- 
hole transformation. 

Now, let us examine how the other terms in the Hamil- 
tonian transform upon the particle-hole transformation 
ip <-» ip^ . For the superconducting pairing potential, we 
find that A — ► A*, so the chirality changes to the op- 
posite. For the impurity potential (11), we find that 
Vi mp — > — Vi m p. On the other hand, the observable Hall 
conductivity a xy should not depend on the choice of op- 
erators ip and ip', which are integrated out. So, we obtain 



Cz:y(A, Vimp) = CT xy (A*, — Vi mp ) 



and, combining with Eq. (18), 



Cxy(A, Vimp) = -&xy(A, —Vimp). 



(19) 



(20) 



Eqs. (19) and (20) arc valid only when the approximation 
(16) is employed, and the electron kinetic energy has the 
particle-hole symmetry. 

Eq. (20) shows that the sign of the Hall conductivity 
must change with the sign change of the impurity po- 
tential. This is indeed the case for the skew-scattering 
diagrams shown in Fig. 3, where a xy oc Vj mp . Thus, we 
conclude that the skew-scattering diagrams can be calcu- 
lated using the linearized approximation (16). The sign 
of the anomalous Hall conductivity given by these dia- 
grams changes when repulsive impurities are replaced by 
attractive ones. In that sense, the sign of the anomalous 
Hall effect produced by skew scattering is not a property 
of the Sr 2 Ru04 material as such, but is a property of the 
impurities in this material. This conclusion is contrast 
to the sign of the conventional Hall effect, which is de- 
termined by the sign of carriers, electrons or holes, but 
not by the sign of the impurity potential. 

On the other hand, the diagrams for the Gaussian 
model of disorder, shown in Fig. 4, give <r xy oc V| mp , 
which is incompatible with Eq. (20). It means that 
these diagrams must vanish when calculated using the 
linearized approximation (16). In order to obtain a non- 
zero result for these diagrams, we must go beyond the 
linearized approximation (16) and take into account the 
curvature of the electron dispersion i.e., to take into 
account the electron-hole asymmetry. This is similar to 
the conventional Hall effect, whose sign is determined by 
the particle-hole asymmetry. In this sense, the anoma- 
lous Hall effect given by the diagrams in Fig. 4 reflects 
the properties of the material, because its sign is deter- 
mined by asymmetry of the electron spectrum and does 
not depend on the sign of the impurity potential. 

It should be emphasized, however, that the magni- 
tude of the anomalous Hall conductivity, given by the 
diagrams in Figs. 3 and 4, is still proportional to some 
power of the strength of the impurity scattering. In this 
sense, the magnitudes of the measured anomalous Hall 
conductivity and the Kerr angle do not characterize the 
Sr2Ru04 material as such, but characterize the degree of 
disorder in the samples of this material. This conclusion 



can be verified experimentally by intentionally introduc- 
ing mild dozes of disorder into Sr2Ru04, such that su- 
perconductivity is not significantly suppressed. The the- 
ory predicts that the observed Kerr angle should increase 
substantially after introduction of disorder. 

Concluding this Section, we remark that the diagrams 
shown in Fig. 2 still vanish, even if we take into account 
the particle-hole asymmetry. Their vanishing is a conse- 
quence of angular integration over electron momentum. 
We also mention that the effects of particle-hole asymme- 
try are important to explain certain properties of the su- 
pcrfluid 3 Hc [61], and, thus, may be relevant for Sr2Ru04 
as well. By taking into account the particle-hole asym- 
metry, Yip and Sauls [39] have shown that a chiral p-wave 
superconductor exhibits circular dichroism and birefrin- 
gence, which arise from the collective modes of the order 
parameter. However, the effect is too small to explain the 
experiment [1]. The particle-hole asymmetry was also 
discussed in relation with the sign change of the Hall ef- 
fect at T c in the high-T c superconductors attributed to 
vortex motion [57-60] . 



III. SKEW SCATTERING IN THE 
NON-GAUSSIAN MODEL OF DISORDER 

In this Section, we study the non-Gaussian model of 
disorder. Within this model, the non-zero contributions 
to the anomalous Hall conductivity a xy are given by the 
diagrams in Fig. 3. These diagrams differ from the dia- 
gram in Fig. 2 by the presence of an additional fermion 
loop. The analytical expressions for the diagrams a) and 
b) in Fig. 3 can be written as 



}x V a) (u n ) = n 3 n iU lTr A x T 3 G(m h p 3 )T 3 A y f 3 



QxyH^n) = KsUiulTr A x T S A y T 3 G(zUl p 3 )f 3 



(21) 
(22) 



where the effective vertices and A y are given by Eqs. 
(14) and (15). The trace is taken over the Pauli matrices, 
the frequency wi , and the momentum p 3 in the additional 
fermion loop. The integration over p 3 yields 



P3 



G{w h p 3 )-- 



-N{0) j 



— in7;fo-/V(0)7r 
^f+F+Af * ^? + Ag ' 

(23) 



Here we used the approximation (16) and introduced the 
2D energy density of states N(0) at the Fermi level 



PF 

2ttvf ' 



(24) 



where V is the volume of the system. For a model with 
the parabolic dispersion en = p 2 /2m e and an effective 
electron mass m e , Eq. (24) gives N(0) = m e /27r, but we 
will write our results for a general dispersion. 
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Substituting Eq. (23) into Eqs. (21) and (22) and eval- 
uating the Nambu traces (Q x 3 y a ^ oc Tr{fiT2T3} = 2i and 
Q x f oc Tr{nT 3 T2} = -2i), we find the following result 
for the full response function Q x 3 y — Q*xy^ + Q x 3 y^ [65] : 



e 2 TrN(0)p 4 F A x A y 



Qxy(Un) = K 3 mUo'- 



H(co n ). (25) 



Here the dimensionless function H(uj n ) is 

1 1 



H(u n ) = A T^\ 



V^'t + A v'( n7 i+<"n) 2 +Ag, 



+ UJ n ) 2 + 



A 2 

o , 



(26) 



The Matsubara sum in Eq. (26) can be evaluated using a 
contour of integration in the complex plane enclosing the 
points of non-analytic behavior, as shown in Appendix C. 
The result is 



+ 



oo 



tanh(f£) 


w„+3iA 


2 




^„(a)„ + 2iA ) 


tanh(^r)(ix 


3ia; — 2u; ra 


V* 2 -a 


2 



_Ag + (ia:-w n ) 



,-3iA 



V / w„(w„-2«A ) 



(27) 



3ix+2ii> 11 



A 2 + (m- 



Then, we perform the analytical continuation in Eq. (27) 

iui n — > w + = w + z<$, (28) 

where u> is the real physical frequency and S is an in- 
finitesimal positive number. Thus we obtain the re- 
tarded finite-temperature response function Q x 3 y (u>) and, 
via Eq. (2), the anomalous Hall conductivity a xy ' (uj) per 
layer 



w 



e 

S X yJ 



WA 



H 



(29) 
(30) 



Here the sign function s xy is given by Eq. (4). We restore 
the factor T% in Eqs. (29) and (30), but assume that the 
frequency u) is measured in the energy units. We also 
introduce the parameter W with the dimensionality of 
(energy) 2 , which characterizes the strength and skewness 
of the non-Gaussian disorder potential. The parameter 
W can be positive or negative depending on the sign of 
the impurity potential uo and will be estimated in Sec. V. 
The dimensionless function H (x) describes frequency and 
temperature dependence of the chiral response 



H(x) 



OC 

I 



tanh 



dy- 



IT 



3y + 2x 3y — 2x 
l-(y+x) 2 ~ l-(y-xf 



1 tanh | ^ 

2 \2T 



+ 



y / x(2-x) y/-x(x + 2) 



Generally, the function H (x) takes complex values when 
the variable x = lv + / Aq changes from — oo to oo. Using 
Eqs. (31) and (29), it is easy to check that a xy {to) sat- 
isfies the causality requirement: <r xy (uj) = <r* y (—uj) [56]. 
Eqs. (29) and (31) give a complete answer for the real 
and imaginary parts of the anomalous Hall conductivity 
<Jxy <{w) for arbitrary frequency and temperature. 

Now we focus on the limit T = 0, where the function 
H(x) can be calculated analytically by setting the tanh 
function to unity in Eq. (31) 



H(x) = 



(x — 3) arccos(l — x) (x + 3) arccos(l + x) 



TTyJx{2 — x) 



TTy/-x(2 + x) 



(32) 



Taking into account that the function arccos(x) of a real 
variable x has a real value for |x| < 1 and an imaginary 
value for > 1, we observe that the second term in 
Eq. (32) is always real for x > 0, whereas the first term 
is real for < x < 2 and imaginary for x > 2. The 
appearance of the imaginary part for x > 2 represents 
the threshold of absorption across the superconducting 
gap for photons with the energies uj > 2A . The real and 
imaginary parts of H{x) exhibit a square-root singularity 
at x = 2: H(x) <~ 1/V2 — x. The other asymptotes are 

-aT, x — > 

7T 

(33) 



H(x) 



1057T 



4 In a; 

— i , x — > oo. 

7TX 



Frequency dependence of the real and imaginary parts 
of a x 3 y \uj) at T = 0, obtained from Eqs. (29) and (32), 
is plotted in Fig. 5. We observe that the imaginary part 
of <r xy appears only for ui > 2A , above the threshold of 
photon absorption. In the high-frequency limit u> ^> Ao, 
a x 3 y((j) is given by Eqs. (29) and (33) 



T (3) 
'xy 



WAp, 



S xy h 



7T LU \Ao 



(34) 



where A is the superconducting gap at T = 0. In the 
dc limit uj 



0, a x 3 J has the finite real value 



<j x V{uj = 0) = s 



xy 



1057T 



e 2 W_ 
¥ A 2 ' 



(35) 



At a finite temperature T^0, A (T) in Eqs. (29) and 
(31) should be understood as the temperature-dependent 
superconducting energy gap obtained by solving the ap- 
propriate BCS equation for the chiral p- wave pairing [61]. 
It can be easily verified that, at high frequency w»T c , 

/ON 

the imaginary part of a xy (uj) is much greater than the 
real part at any T < T c , as in Eq. (34) at T = 0. Thus, 
we focus on temperature dependence of the imaginary 
(31) part a x f \u>,T). Taking the limit x > 1 in Eq. (31) and 
observing that the first, integral term gives a negligible 
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s 



-10 L 

FIG. 5: (Color online) Frequency dependence of the real 
(dashed line) and imaginary (solid line) parts of the anoma- 
lous Hall conductivity a^J (lu) at T = given by Eqs. (29) 
and (32). 




T 



0.2 0.4 0.6 0.8 1.0 T c 

FIG. 6: (Color online) Temperature dependence of the imag- 
inary part of the anomalous Hall conductivity (Jx^\lu,T) for 
u > T c given by Eq. (36) with A (T) from the BCS theory 
[67]. 



contribution to the imaginary part, we find from Eq. (29) 



t^ 3) (uj,T) 



_ Aq(T) 

(w,0) A o(0) 



tanh 



Ao(T) 
IT 



cj>T c . (36) 



The plot of Eq. (36) is shown in Fig. 6 using A (T) cal- 
culated from the BCS equation [66, 67]. For T close to 
T c , where Aq(T) <C T, the tanh function in Eq. (36) can 

be replaced by its argument, so a xy becomes propor- 
tional to the square of the superconducting order param- 
eter Aq(T). In the Ginzburg-Landau theory, Aq(T) oc 
(T c — T), so Eq. (36) gives a linear temperature depen- 
dence for a x 3 J (and, thus, 6k) near T c . The factor A 2 , 
originates from the product A x A y in Eq. (25), because 
both A x and A y are necessary for a non-zero anomalous 

Hall effect, as emphasized in Ref. [31]. However, a x 3 J at 
T = in Eq. (34) is proportional to the first power of Ao, 
somewhat reminiscent to what was proposed in Ref. [1]. 
Curiously, the temperature dependence in Eq. (36) is 
the same as for the critical current in Josephson junc- 
tions [68]. 

For high frequency lu 3> T c and T close to T c , we obtain 
the following expression for the real and imaginary parts 
ofo-$$(u,T) from Eqs. (29) and (31) 



r(3) 



e 2 WAp(T) 
s xy% 



AT r 



■In 



(37) 



The imaginary part a 1 ' oc 1 /lu 3 in Eq. (37) looks similar 
to the real part of a xy calculated by Goryo in Eq. (5) of 
Ref. [48] for the same model in the same asymptotic limit. 
Given the causality requirement a xy {io) = a* y (—Lu) [56], 
it appears that the factor of i = was overlooked 

in the calculation of Ref. [48] , so the real and imaginary 
parts of a xy were interchanged. 



IV. MODEL OF THE GAUSSIAN DISORDER 
WITH PARTICLE-HOLE ASYMMETRY 

In this section, we calculate the anomalous Hall con- 
ductivity for the Gaussian model of disorder by taking 
into account the effects of particle-hole asymmetry. 

The lowest-order non-vanishing diagrams in this case 
are shown in Fig. 4. Fig. 4b actually represents four 
similar diagrams, with self-energy corrections to the up- 
per and lower Green's functions. These diagrams can be 
considered as the lowest-order terms of a more general 
diagrammatic series shown in Fig. 7, which is often con- 
sidered in calculations of transport properties of metals 
with disorder. We find it more practical to calculate the 
general series of diagrams in Fig. 7 and then take the limit 
of low concentration of impurities, which corresponds to 
the diagrams in Fig. 4. In this way, we can be sure that 
all contributions of the same order are taken into account, 
and mutual signs of the diagrams are correct. 



A. Self-consistent non-crossing approximation 

We begin with calculation of the Green's function 
Qp{wi) dressed due to impurity scattering. (Starting 
from this Section, we write the momentum p in the ar- 
gument of a Green's function as a subscript in order to 
shorten notation in long mathematical equations.) The 
dressed Green's function is obtained from Dyson's equa- 
tion shown in Fig. 7a, where we make the standard non- 
crossing approximation, assuming that Epr 1 and ne- 
glecting diagrams with intersecting impurity lines. The 
self-energy ~S(wi) due to impurities is 

= niul^Gpi-wi) = iwi[l - T?i(tz7i)]f - mi^i)^- 

(38) 
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a) 



Go 



+ 



X 



Go So Q 



6) 



as we will see later, this divergence cancels out in the final 
results, and the outcome can be expressed in terms of 
N'(0), which is a measure of the particle- hole asymmetry. 

Using S(o7;) from Eq. (38) to solve Dyson's equation 



we obtain the dressed Green's function 



(43) 



c) 




FIG. 7: Diagrams for the disorder- dressed Green's function 
(a), the vertex function (b), and the response function (c). 



One can notice that the terms with fi and f 2 vanish after 
angular integration of A(p) over p in Eq. (38). Here we 
introduced the functions r]i(wi) and 772(07;) defined as 



7/1 (to; ) = 1 



1 



AT 



1 f 

m{™{) = — / 

T7T J_ 



100 N(Z)dt 



N(0) rf + At+i 



2 ' 



(39) 

(40) 
(41) 



where r is the quasiparticle lifetime due to scattering 
on impurities. One can notice that, within the Gaus- 
sian model of disorder, the concentration of impurities 
rii and the strength of the impurity potential Uq can be 
completely absorbed into the definition of the quasipar- 
ticle lifetime r in Eq. (41). This is not the case for the 
non-Gaussian model of disorder discussed in Sec. III. 

The function r]i(wi) in Eq. (39) is well known in the 
theory of impurity scattering. Moreover, the first term 
with fo in Eq. (38) is, essentially, the same as in Eq. (23). 
In contrast, the function rj 2 (rui) in Eq. (40) is less com- 
mon and deserves more discussion, because it will play 
crucial role in this Section. The function iV(£) in Eq. (40) 
is the energy density of states as a function of the electron 
energy £. If we make the approximation iV(£) = N(0), 
i.e., take N(£) to be a constant, then the integral in 
Eq. (40) vanishes, because the integrand is odd in £. 
However, if we take into account that N(£) is not a 
constant, i.e., consider the particle-hole asymmetry, then 
r) 2 (m{) is non-zero. Let us expand N(£) near the Fermi 
energy 



N(0~N(0) + N'(0)^ 



(42) 



where iV'(0) = dN/d£ is taken at £ = 0. When Eq. (42) is 
substituted into Eq. (40), the integral diverges. However, 



iwir\ X Tv + (Cp -ri2)h +p x A x f 1 -p y A y T 2 



(44) 



In order to calculate the two-particle response func- 
tion sclf-consistently, it is necessary to include vertex 
corrections due to impurity lines. They transform the 
bare vertex Jj(p) — — eVj(p)ro of interaction with the 
electromagnetic field in Eq. (6) into the dressed vertex 
Tj (p) . The vertex Tj (p) is obtained by solving the Bcthc- 
Salpeter equation illustrated in Fig. 7b 

Fj(p) = 7j(p) +n.«^T35 J) '(n'i)fj(p')5p'(ni+w„)T3. 
P > 

(45) 



For elastic impurity scattering, the integral equation (45) 
can be solved analytically by summing the geometric se- 
ries as shown in Appendix D. The solution for T y , which 
contributes to the anomalous Hall response [see Eq. (52)], 
is given by Eq. (D4) 



T y (p) = ly{p) + Ay P pVF (-iain + a 2 f 2 )L(m u u n ) 



I-61 



2ttt 
i-Tb 



ib 2 



(l-6i) 2 + 6l u (l-foi) 2 + &| 



T3 



(46) 



Here the functions a\ and a 2 are defined as 



a-i{wi,u} n ) = rfrim) ~ mfai + u n ), (47) 
a 2 (zui,uj n ) = iwir)i(zui) + w n ) -q^wi +u> n ), (48) 

and the function L(wi,uj n ) is 



L(zUi,LO n ) 



/uin 
-UJn 



N(Z)d£ 



_ ; . N(0) D{wi,S)D{wi+W n ,Z) 

( 1 1 \ 



D(m h Z) = rfvKm) + K - rni^i)] 2 + A 2 . 



(49) 
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The functions b\ and hi are given by the integrals 



2 f"D 



N(0) 

[£_-r 1 2(u7 l )][^-r]2(TAj+)}+miw+T] 1 (nJi)i ll (w+) 



(50) 



D(o7,,0^i + .0 



6 2 {m 



1 r 



CT ; + 7?i(CT ; + )[g-T?2 (OT;)]-OT;T?i(OT;)[g-T? 2 (^ + )] 
D(o7,,0^i + »0 



(51) 



where tu ; + = n7; + u) n . 

Using Eqs. (44) and (46), one can now calculate the 
dressed response function Q x 4 y shown in Fig. 7c [65] 

fi$K) = Tr [7x(p) r„(p) S P (^ + w„)]. (52) 

This expression can be simplified by integrating over the 
angle of p first. After substituting the dressed vertex 
T y (p), the contribution from the first term in Eq. (46) 
drops out after taking the trace over the Nambu space, 
as discussed in Sec. II B. Thus, we are left only with the 
second, momentum-independent term in Eq. (46), and 
the integral over p in Eq. (52) reduces to 



^ Gp {-CD l + OJn ) % (p) Gp{wi) 



47T 



(a 2 T 1 - mif 2 ) L(m h uj n ). 



(53) 



[One can notice that, in the zeroth order in disorder (771 = 
1 and r]2 = 0), Eq. (53) reduces to Eq. (14).] Substituting 
Eq. (53) into Eq. (52), i.e., multiplying Eq. (53) by the 
second term in Eq. (46), and taking the trace, we obtain 
the response function 

(4) e 2 A x A yP 3 F v F ^ 2 



-4iaia 2 (l - 61) - 2 (a? + a 2 .)^ 



(l- &1 )2 + 6 2 



In contrast to the clean case discussed in Sec. II B, the re- 
sponse function (52) does not vanish identically, because 
the second term in Eq. (46) is momentum-independent 
as a consequence of randomization of electron momen- 
tum due to scattering on impurities. However, if the 
model has the particle-hole symmetry, then the functions 
a\ (47) and 6 2 (51) vanish, and, thus, the anomalous Hall 
response (54) is zero to any order in impurity scattering. 
To obtain a non-zero result, we need to take into account 
the particle- hole asymmetry explicitly, as discussed in the 
next Subsection. 



Perturbative expansion in the strength of 
disorder 



As mentioned in Sec. II C, Sr 2 Ru04 is a stoichiomet- 
ric crystalline material with a very low concentration of 
impurities, and the existence of p-wave pairing requires 
that Ao S> 1/t. Thus, we can use the strength of disor- 
der as a small parameter and expand the anomalous Hall 
response function (54) to the lowest non-vanishing order 
in 1/t. The non- vanishing term appears in the order 
(1/t) 2 , i.e., in the fourth order in Vj mp : 

Wl 

x {2wi + uj n ) [rfr{m) - mi^i + w n )]. 

Notice that one power of 1/t appears in the prefactor in 
Eq. (55), and another power comes from the function 772 
defined in Eq. (40). Eq. (55) precisely corresponds to the 
diagrams shown in Fig. 4. 

Using Eqs. (40) and (42), the difference 772 (a?;) — 
filial + w„) in Eq. (55) can be written as 



1 



N(0) 



(56) 



1 N'(0) 



t N(0) 

2 jV'(0)a>„(2n7, +^») 

TIT N(0) U1 D 



(57) 



When Eq. (42) is substituted into Eq. (56), the constant 
term N(Q) in the density of states cancels out, and the 
term proportional to N'(0) produces a convergent inte- 
gral over £. Thus, the anomalous Hall response (55) be- 
comes directly proportional to the particle-hole asymme- 
try parameter N'(0) and would vanish if N'(0) = 0. 

For generality, we consider the finite limits of inte- 
gration ±uod in Eq. (56). The first term in Eq. (57) 
represents the value of the integral when the limits are 
taken to ±00, whereas the second term is a correction 
due to the finite limits of integration, assuming that 
uod 3> A ,uj n . As we will see below, one of these two 
terms gives the dominant contribution to the imaginary, 
and the other one to the real part of the ac Hall con- 
ductivity in the high-frequency limit. We first evaluate 
the cutoff- independent contribution to a xy {uo) given by 
the first term in Eq. (57) and then the cutoff-dependent 
contribution from the second term. 



1. Cutoff-independent contribution to the anomalous Hall 
conductivity 

Substituting the first term in Eq. (57) into Eq. (55) 
(or taking the limit lojj — > 00), we obtain the cutoff- 
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independent contribution Q^y to the anomalous Hall 
response function 



3e 2 i/A 



)xy {Un) - Sxy ^2^-2^2 



K(cj n ), (58) 



where the dimensionless parameter v is a measure of the 
particle-hole asymmetry 



_ p F v F N'(0) 
V ~ 2N(0) ' 

and the dimensionless function K(ui n ) is 
TA ( 



(59) 



1 



^ 2nj; + w„ ^ ^/Wf+Ag v / (ro i +w„) 2 + A 2 / 

(60) 



x ,/(w n + w/) 2 + A 2 - Jwf + A 2 



The Matsubara sum in Eq. (60) can be evaluated using 
complex analysis as shown in Appendix E. The function 
K(u) n ) can be written as the sum of three terms 

K(uj n ) = Kx(u n ) + K 2a (u n ) + K 2b (uj n ), (61) 

where the functions Ki(u> n ), K2 a {w n ), and K2b(w n ) are 
given by Eqs. (E4), (E6), and (E7). After the analytical 
continuation (28), we obtain 



dx 



-H2u 



27T y/ X 2 - 1 (d> 2 - 4.T 2 



■ tanh 



xA 
2T 



K 2 aM = ~tanh 



K 2b {uj) 



f°° dx i-z 

2i — V x 2 - 1 tanh 

Jl 27T 
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xA Q 



2T 



(62) 

(63) 
(64) 



1 



.(2a?+w+)(l - [x+Q+] 2 ) (2x-Q+)(l - [x-u+] 2 ) 



where lu + = (uj + i5) /Aq. 

Using Eqs. (58) and (2), we find the cutoff-independent 
contribution a^f* to the anomalous Hall conductivity 



*^W^(!sD' (65) 



T (4a), 



^t a) (o) 

10 r 



.. — r" 



FIG. 8: (Color online) Real (dashed line) and imaginary (solid 
line) parts of the cutoff-independent contribution a^'(ui) to 
the anomalous Hall conductivity at T — given by Eqs. (65) 
and (66). 



The plots of the real and imaginary parts of <j X y vs. u> 
at T = 0, obtained from Eqs. (65) and (66), are shown in 
Fig. 8. At y — 2, the function K(y) exhibits a square-root 
singularity. The asymptotes of K(y) are 



-iK(y) 



15tt 



y^O 



(67) 



y^oo. 



So, at high frequencies u> Aq, a X y becomes 
- 2 3z/A 



and the dc limit is 



i H hi . . 

7T CJ \ Aq 



4t a) (^ = o) 



e z/ 



S 57T 3 V tAq 



1 



(68) 



(69) 



The electron-hole asymmetry factor v given by Eq. (59) 
is, generally, of the order of unity, whereas the factor 
rAo 3> 1 in Eq. (69) is large, as required for the existence 
of p-wave superconductivity. Thus, the anomalous dc 
Hall conductivity (69) per layer is always much smaller 
than the quantum of conductance e 2 /h. 



where the dimensionless function K (61) is given by the 
sum of Eqs. (62), (63), and (64). 

At T = 0, the function K(y) becomes 



6i arcsin(y/2) 



y-2J ir y/4-y 2 



(66) 



+ i 



dx 
n 



(2x+y)(l - [x+y] 2 ) (2x-y)(l - [x-y] 2 ) 



2. Cutoff- dependent contribution to the anomalous Hall 
conductivity 



Now we consider the contribution Q X y^ to the anoma- 
lous Hall response function (55) from the second term in 
Eq. (57), which depends on the energy cutoff lod- 



6e 2 vA 



-J{Un), 



(70) 
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where the dimensionless function J(u> n ) is 



i \ V m i ~ 



1 



a 2 ^{w l +uj n y+/^ 



(71) 



The Matsubara sum in Eq. (71) is calculated in Ap- 
pendix F. The result for J(uj n ) is 



J(w n ) = tanh 



Ao 
2T 



dx 



Vx 1 



A 2 



tanh (— ^ 



x 
2TJ 



Ar 



A 



(72) 



\^{ix + uj n ) 2 + A 2 yj{ix - uj n ) 2 + A 2 ^ 

Performing the analytical continuation (28) in Eqs. (72), 

(70), and (2), we obtain the cutoff-dependent term cr^f 1 
in the anomalous Hall conductivity: 



6^A 



xy h t: 3, t 2 uj d uj 2 " ^A 
Here the dimensionless function J(jj) is 

tanh 



(73) 



J(y) = tanh('^-- 



2T 
1 



dx 



V x 2 



1 



xA 
2T 



1 



y/1 - (x - y) 2 ^l-(x + y y 



(74) 



At T = 0, the explicit analytical expression for the 
function J{y) can be obtained from Eq. (74): 



oo 



1 



1 



Va; 2 -! Wl-(x-y) 2 v/l-(a: + y) 2 



-1 - 



TTV^ 



F\ iArcsinh 



F zArcsinh< \ / — 



y 



(75) 



The function F(x; k) is the incomplete elliptic integral of 
the first kind. The plots of the real and imaginary parts 
of critf' 1 vs. ui at T = 0, obtained from Eqs. (73) and (75), 
are shown in Fig. 9. At y = 2, the imaginary part of J(y) 
has a discontinuity (jump), whereas the real part of J(y) 
has a logarithmic divergence. The function J(y) has the 
following asymptotes: 



J(y) = 



16 



V 



y^Q 



4i 1 
ny 



(76) 



oo 



So, at high frequencies a; 3> Aq, Oa^f becomes 
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FIG. 9: (Color online) Real (dashed line) and imaginary (solid 
line) parts of the cutoff-dependent contribution critf' (w) to the 



anomalous Hall conductivity at T 
(75). 



and the dc limit is 



given by Eqs. (73) and 



4 4 JV = o) 



6v 



h 16tt 3 t 2 ojdA 



(78) 



3. Comparison of the cutoff- dependent and -independent 
contributions to the anomalous Hall conductivity 

Combining Eqs. (65) and (73), we obtain the total 
anomalous Hall conductivity o~ X y within the Gaussian 
disorder model with particle-hole asymmetry 



a W = CT (4a) + (46) 
xy xy 1 xy 



(79) 



At high frequencies uj 3> Ao, the asymptotic expres- 
sions for the real and imaginary parts of the anomalous 
Hall conductivity (79) are 



(Uj) — S X y- 



3vA 



3^A 



xy 



ujduj 

1 

i3 



6Ao 
• 4 V Ao 



to ' 
8A 



uj-' 



TT^UJrjUJ 



3 1 A 



(80) 



(81) 



Eq. (81) shows that the dominant contribution to the 
imaginary part o'^y^ is given by the cutoff-independent 
term. In contrast, Eq. (80) shows that, in the case of 
Ao<w< ll>d <C w 2 /[Ao ln(tj/Ao)], the dominant con- 
tribution to the real part o~' X y is given by the cutoff- 
dependent term. This is why it is important to perform 
calculations keeping the cutoff u>jj large, but finite. 

On the other hand, the low-frequency anomalous Hall 
conductivity is dominated by the cutoff-independent 
term: aiffl ^> a X y from Eqs. (69) and (78) at uj = 0. 
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EXPERIMENTAL IMPLICATIONS 



the non-Gaussian model 



A. Estimates of the anomalous Hall conductivity 

In this Subsection, we discuss experimental implica- 
tions of our results (29) and (79) for the anomalous Hall 
conductivity. First, we estimate the magnitude of the 
anomalous dc Hall conductivity a xy (to = 0). 

Within the Gaussian disorder model, the dc Hall con- 
ductivity is given by Eq. (69). [As discussed in the 
preceding Section, Eq. (78) gives a negligible contribu- 
tion compared with Eq. (69).] For a crude estimate, we 
use the following values of the parameters in Eq. (69): 
Ao = 0.23 meV, as estimated in Sec. IB from the value 
of T c , 1/r w 7x 1(T 5 eV taken from Ref. [1], and v ~ 10. 
The particle-hole asymmetry parameter v (59), in prin- 
ciple, can be deduced more accurately from the band- 
structure calculations [69] for Sr 2 Ru04. Thus we obtain 
an estimate for a x 4 J from Eq. (69) 



7 (3) 

'xy 



0),18-. 



(84) 



(w = 0)~HT 2 V- 



(82) 



Because tAo 3> 1, the dc Hall conductivity (69) within 
the Gaussian disorder model is much smaller than the 
conductance quantum e 2 /h. 

Now we consider the anomalous dc Hall conductivity 
a xy (u) = 0) originating from the non-Gaussian disorder 

model. Eq. (35) for a xy (uj = 0) contains the parameter 
W (30) proportional to the combination n^riiV^, which is 
different from the combination n^u 2 , appearing in Eq. (41) 
for r. Thus, the knowledge of the impurity scattering 
time t alone is not sufficient to calculate W . In addition 
to r, one needs to know the concentration of impurities rij 
and the skewness parameter K3. Goryo [48] estimated the 
typical distance between impurities as k = 1000—5000 A. 
For the estimates in our paper, we use the value k — 
1000 A, which corresponds to rij ~ 10 14 m~ 2 [70]. Using 
Eq. (24) for N(0) and Eq. (41) for 1/r, we rewrite the 
expression for W in the following form 



\W\ 



«3 



/ p 3 F hv F 

32n,T 3 



(6 meV) 2 



(83) 



For the numerical estimate given in Eq. (83), we used 
the values p F = 0.75 A -1 and v F = 5.5 x 10 4 m/s 
from Ref. [5] (giving Ep — ppvp/2 = 0.14 eV) and 
the values of 1/r and nj quoted above. For an esti- 
mate of the skewness parameter, we took the upper limit 
K3 = 1, not having information about the actual nature 
of impurities. This value corresponds, for example, to 
randomly distributed impurities generating short-range, 
delta-function potentials of the equal strength uq, as dis- 
cussed after Eq. (17). The same assumption was utilized 
in Ref. [48]. 

Substituting Eq. (83) into Eq. (35), we obtain an esti- 
mate for the magnitude of the dc Hall conductivity within 



Eq. (84) gives a value much greater than the conductance 
quantum e 2 /h, even though rA 3> 1. Thus, for the pa- 
rameters given above, we find that the contribution from 
the skew-scattering diagrams to a xy (uj = 0) is dominant. 
This is because the skew-scattering processes correspond 
to the lower-order diagrams in impurity concentration. 
Finally, we emphasize that our calculation is done for 
a single domain of the p x ± ip y superconductor, while 
a realistic macroscopic sample of Sr 2 Ru04 should con- 
sists of multiple domains with opposite chiralitics. The 
dc Hall effect contributions of the opposite signs from 
different domains would cancel out and make an experi- 
mental verification of Eq. (84) difficult. 

We now calculate the anomalous ac Hall conductivity, 
given by Eqs. (29), (65), (73) and (79), at the optical fre- 
quency ui — 0.8 eV utilized in the experiment [1]. Some 
of these equations contain the cutoff frequency up, , which 
depends on the microscopic nature of the pairing inter- 
action and is not known very well. In our calculations, 
we assumed that u>d > Using the optical frequency ui 
as the lower bound for the cutoff frequency wd, wc find 



erg) (to = 0.8 eV) - (10~ 10 + 



(85) 



41? (w = 0-8 eV) ~ (10- 12 + lO- 1 ^) —. (86) 

Comparing Eqs. (85) and (86), we observe that the non- 
Gaussian model gives a much greater contribution than 
the Gaussian model to both real and imaginary parts of 
<j xy (oj — 0.8eV). Thus, in the rest of this Section, we will 
use Eq. (85) for the non-Gaussian model to estimate the 
Kerr angle. However, one should keep in mind that the 
dominance of (Jxy IS the consequence of the high estimate 
for the parameter W in Eq. (83). This estimate depends 
on the parameters K3, rn, and uq, for which there are 
no direct measurements. Similarly, a x 4 J for the Gaussian 
model depends on the parameters v and ojp>, for which 
there are no direct measurements cither. The parame- 
ter r, utilized for both models, is only indirectly inferred 
from the dc conductivity in Ref. [1]. Given these great 
uncertainties, our numerical estimates should be consid- 
ered as only tentative. 



B. Estimates of the polar Kerr angle 

The anomalous Hall conductivity a xy contributes to 
the dielectric permeability tensor *e, which character- 
izes propagation of electromagnetic waves in the medium. 
The dielectric permeability tensor *e is related to the con- 
ductivity tensor a as follows 



e = £00 H o. 



(87) 
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Here is the background dielectric tensor, which orig- 
inates from polarizability of the other, non-conduction 
bands in the material. The diagonal components of the 
conductivity tensor are assumed to have the Drude-like 
form: 



(Tjj = - 



(88) 



where uij and jj are the plasma frequency and the quasi- 
particle scattering rate along j-th axes, respectively. Due 
to the square symmetry of Sr 2 Ru04 in the afe-plane, we 
have u x =u y = us ab , "f x = y y = j a b, and e xx = e yy = e ab , 
so 



£ a6 (w) = £ c 



J ab 



(89) 



It should be emphasized that we distinguish between 1 /r, 
the scattering rate on impurities given by Eq. (41), and 
7, the quasiparticle scattering rate at optical frequencies, 
which may be dominated by other scattering processes. 

Now let us consider a polarized electromagnetic plane 
wave incident in the z direction normal to the (x, y) sur- 
face of a Q2D chiral superconductor, as relevant for the 
experiment [1] . Propagation of the electromagnetic wave 
in the medium is described by Maxwell's equation: 



E = 0, 



(90) 



where E is polarized in the (x, y) plane. The anoma- 
lous Hall conductivity a xy produces antisymmetric off- 
diagonal matrix elements e xy = — e yx in the tensor *e 
via Eq. (87). Thus, propagation of the wave inside of 
the superconductor is described by two circularly polar- 
ized cigenmodes with different refraction indices n + and 
n_. Given that the off-diagonal elements are typically 
very small \e xy \ <C \e xx \, one can expand n + and n_ to 
the first order in e xy oc a xy . Then, using the boundary 
conditions for the electric field at the interface between 
vacuum and the material, one finds the reflection coeffi- 
cient \r\ and the polar Kerr angle 9 k [24, 36] 



l?1 |n+l|' 

4-7T 



l (™ 2 - 1) yj£ab{u) [e a &M - 1] 



(91) 
(92) 
(93) 



Here n(u) = \/e a b(oj) is the complex refraction coeffi- 
cient. In Eq. (92), a xy (uj) is the 2D Hall conductivity 
per one layer calculated in the previous Sections, and 
the interlayer distance d — 6.8 A [5] converts it into the 
three-dimensional, bulk conductivity implied in Eq. (87). 
As Eq. (92) shows, the magnitude of the Kerr angle 9k 
depends not only on the anomalous Hall conductivity 
a xy (ui), but also on the complex refraction coefficient 
n(u>). The frequency dependence of the refraction coeffi- 
cient is beyond the scope of the present paper, in which 



we focus on the TRSB effects. Therefore, we estimate 
the parameters determining n(u>) from the experimental 
data [71], where the out-of-plane optical conductivity of 
Sr 2 Ru04 was studied. Because many parameters (e.g., 
the quasiparticle scattering rate 7) for the diagonal com- 
ponent a xx (u) of the conductivity tensor are difficult to 
estimate, it is desirable to measure n(w) experimentally 
for the same samples where the Kerr effect is measured. 

According to Ref. [71], the dielectric constant and the 
plasma frequency in the a&-plane are £00 = 10 and w & = 
2.9 eV [72]. Substituting these values into Eq. (89), we 
find the frequency of the plasma edge lo p — uj a b/ '\/ £ '0° = 
0.9 eV. So, the measurement frequency u> = 0.8 eV in 
the experiment [1] appears to be below the plasma edge 
frequency lu p . Another important parameter here is the 
quasiparticle scattering rate 7 a &. There is experimental 
evidence that 7^ is quite large, of the order a fraction of 
eV [71]. To illustrate the importance of this parameter 
for the estimate of 9k, below we discuss two limits: u » 
j ab and uj ~ 7 a6 . 

The case of w > "f a b was previously considered by 
Goryo [48], who obtained the value for the Kerr angle 
to be about 30 nrad. However, as we showed in Sec III, 
the analytical structure of the Hall conductivity obtained 
by Goryo is incorrect, which leads to a much smaller 
estimate for the Kerr angle as we discuss below. Indeed, 
consider the expression (92) for the Kerr angle in the 
following form: 



>k = 



47T 

ud 



[a>K>)+a'(u,K>)] 



(94) 



9 K - 



where a'(ui) and a"(u) are the real and imaginary parts 
of a(uj). For the frequency lu < uj p below the plasma 
edge, Eq. (89) gives e ab (to) < 0, so the index of re- 
fraction ti(lu) = y^Sabiw) and a(uj) in Eq. (93) are 
imaginary. Therefore, the main contribution to 9k in 
Eq. (94) comes from real part of a xy . Using the estimate 
= 0.8 eV) - l0- w e 2 /h from Eq. (85), wc find 
O.f nrad. This estimate is much smaller than the 
one obtained by Goryo, because the real part of a xy is 
three orders of magnitude smaller than imaginary part 
that was used in his estimate. 

We now consider the other limit ut ~ 7^. The experi- 
mental measurements of "f a b in the normal state at T — 
9 K by Katsufuji et al. [71] indicate that the quasipar- 
ticle scattering rate at the optical frequency u> ~ 0.8 eV 
is 7a& ~ 0.4 eV. Assuming that the magnitude of the 
quasiparticle scattering rate is the same in the supercon- 
ducting state at T < T c , the real and imaginary parts 
of a(u>) become of the same order as shown in Fig. 10. 
Then, according to Eq. (94), both real and imaginary 
parts of the Hall conductivity (85) give comparable con- 
tributions to the Kerr angle. Using Eq. (89) with the 
value 7a& w 0.4 eV quoted above and Eqs. (93), (94), 
and (85), we estimate the Kerr angle in this case as 



9k ~ 43 nrad. 



(95) 
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FIG. 10: (Color online) Real and imaginary parts of the 
dimensionless function a.{u) given by Eqs. (93) and (89). 
The solid (blue) and dotted (black) lines represent a'(u)) for 
7 = 0.4 and 0.1 eV, respectively. The dashed (red) and dash- 
dot (brown) lines represent a"(ui) for 7 = 0.4 and 0.1 eV, 
respectively. The parameters u> a t = 2.9 eV and £00 = 10 are 
taken from Ref. [71]. 



This estimate is of the same order of magnitude as the 
experimentally observed 9°^ s» 65 nrad [1]. 

In the process of calculating the estimates, we made 
numerous approximations and assumptions. The values 
of many parameters that appear in the theory are un- 
known, and more measurements are needed to determine 
them. Nevertheless, the final estimate (95) is encourag- 
ing and indicates that impurity scattering may provide a 
viable explanation for the polar Kerr effect in S^RuO^ 

When a' and a" in Eq. (94) are of the same order, 
&k is dominated by a'J. y , which is much greater than 
a' xy . Thus, 0k{T) cx a xy (T), so the temperature depen- 
dence of cr xy (T) determines the temperature dependence 
of 9k{T), assuming that a at high frequencies does not 
depend significantly on temperature for T < T c . Then, 
Eq. (36) and the plot in Fig. 6 give the temperature 
dependence of the normalized Kerr angle 9k(T)/9k(0). 
This theoretical result can be directly compared with the 
temperature dependence of 9k{T) experimentally mea- 
sured in Ref. [1]. 



VI. CONCLUSIONS 

In this paper, we calculate the anomalous (sponta- 
neous) ac Hall conductivity a xy (u>) for a chiral p x + ip y 
superconductor, such as Sr2Ru04, in the long wavelength 
limit q — ► in the absence of an external magnetic field. 
We show that the anomalous Hall conductivity vanishes 
for a translationally-invariant system, and a non-zero re- 
sult requires presence of impurities. Non-zero contribu- 
tions to the anomalous Hall conductivity appear in the 
higher order (above the second order) of the perturbation 
theory in impurity scattering. We consider two mod- 
els of disorder potential: Gaussian and non-Gaussian. 
The Gaussian model is characterized by the second mo- 
ment of the random impurity potential, whereas the non- 
Gaussian model has a non-zero third moment. For the 
Gaussian disorder model, we present a symmetry argu- 
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ment and demonstrate by direct calculations that a non- 
zero anomalous Hall conductivity requires particle-hole 
asymmetry of the electron spectrum. Thus, the magni- 
tude and the sign of a xy (u) depend on the band-structure 
parameter N'(0), which is the derivative of the normal 
density of states at the Fermi level. On the other hand, 
the anomalous Hall conductivity in the non-Gaussian 
model is not proportional to the particle-hole asymme- 
try parameter, and its sign is determined by the sign 
of the impurity potential. Therefore, in principle, these 
two models can be distinguished experimentally by in- 
troducing positively and negatively charged impurities. 
However, in practice, such a discrimination may be diffi- 
cult, because the sign of the anomalous Hall conductivity 
also depends on the sign of the chirality of the supercon- 
ducting p x ± ip y order parameter. 

By calculating the lowest-order non-vanishing Feyn- 
man diagrams shown in Figs. 3 and 4, we obtain closed- 
form expressions for the frequency and temperature de- 
pendences of the anomalous Hall conductivities originat- 
ing from these two models. As a function of frequency, 
the calculated Hall conductivities a xy (u) have a finite 
real value at oj — and exhibit singularities at the thresh- 
old of photon absorption uj = 2Ao. At high frequencies 
lj Ao, the dominant contribution to the anomalous 
Hall conductivity comes from the imaginary part a' xy (io) 
which decays as 1/w 3 , see Figs. 5, 8, and 9. As a function 
of temperature T, the high-frequency a''(u>, T) increases 
linearly with the decrease of temperature near T c and 
saturates at T — ► 0, as shown in Fig. 6. 

In order to estimate whether the Gaussian or non- 
Gaussian term gives a dominant contribution to the 
anomalous Hall conductivity for S^RuGi, it is neces- 
sary to know the parameters characterizing the electron 
spectrum and the strength of the disorder. The latter 
are poorly known. Nevertheless, by making a number 
of assumptions and approximations, we obtain numeri- 
cal estimates, which show that the non-Gaussian term 
dominates over the Gaussian one both at high and low 
frequencies. 

The polar Kerr angle 9k is proportional to the ac Hall 
conductivity a xy {u>). However, the proportionality rela- 
tion also involves the complex refraction coefficient n(u>), 
which is poorly known for Sr2Ru04. Using a Drudc 
model with the experimentally estimated parameters to 
obtain n(uS) and the anomalous Hall conductivity from 
our calculations, we estimate the magnitude of the polar 
Kerr angle at u = 0.8 eV as 9k ~ 43 nrad, which is of 
the same order as the experimentally observed value of 
65 nrad [1]. Despite numerous assumptions and approx- 
imations used, this result is encouraging and indicates 
that impurity scattering models may provide a viable 
explanation for the polar Kerr effect in Sr2Ru04. How- 
ever, the main conclusion of our paper is not a particular 
value of 9k, but the qualitative understanding that the 
magnitude of the anomalous Hall conductivity and the 
Kerr angle are directly proportional to some power of the 
impurity concentration and should strongly vary among 
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different samples. Systematic measurements of the Kerr 
angle as a function of concentration of the intentionally 
introduced impurities would be very desirable. The study 
should be performed for sufficiently low concentration in 
the range A » 1/r, where impurities do not affect the 
superconducting order significantly, but do affect the po- 
lar Kerr effect. 
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APPENDIX A: THE PEIERLS-ONSAGER 
SUBSTITUTION FOR SUPERCONDUCTORS 

The BCS mean-field action for a superconducting sys- 
tem is given in Eq. (5). The corresponding Hamiltonian 
H(p) has the form of a 2 x 2 Nambu matrix 

h ^-{a*(p) -z{-p))> (A1) 

acting on the spinor [4>(p), ^ {— p)\- We omitted the 
spin indices of the fermion operators, because they are 
not essential for the consideration here. The notation in 
Eq. (Al) is the same as in Sec. II A. For simplicity, we 
will assume that £(— p) = which is the case for most 
centrosymmetric materials, including Sr 2 RuC>4. 

When the system is subject to an electromagnetic 
field, the question arises how the vector potential A 
should be introduced into the Hamiltonian (Al). For 
non-superconducting systems, it is introduced via the 
Peierls-Onsager substitution p — > p—eA, and it is tempt- 
ing to make the same substitution in Eq. (Al). How- 
ever, one should be careful and re-examine how this 
substitution originates from the requirement of gauge 
invariance. The two diagonal terms in Eq. (Al) give 
the following contributions J d 3 r ijp (r) £(— iV) ip(r) and 
/ d 3 r ip(r) £(— iV) ip^(r) to the Hamiltonian of the sys- 
tem in the real-space representation. When we make a 
gauge transformation of the fermion operators ip(r) — > 
ip{r) e tv ( r \ the gradients of the phase ip need to be com- 
pensated by the gauge transformation of the vector po- 
tential A. This leads to the substitution p — > p — eA 
and p — > p + eA in the upper and lower diagonal term in 
Eq. (Al). The substitutions are different because of the 
different order of ip and ift in these terms. This is very 
well known [43]. 



For unconventional superconductors, where the pairing 
potential A(p) explicitly depends on the electron momen- 
tum p, the question arises whether the Peierls-Onsager 
substitution should be made in A(p). However, it is not 
clear whether the rule p — > p— eA or p — > p + eA should 
be used in the off-diagonal term in Eq. (Al). Actually, 
the Peierls-Onsager substitution is not needed in the off- 
diagonal term, because, when properly written, it is al- 
ready gauge-invariant without introduction of the vector 
potential A. When the order parameter varies in space, 
the pairing potential (3) should be written as the sym- 
metrized combination, i.e., the anticommutator, of the 
two-component order parameter \t' = (A x ,iA y ) and the 
momentum operator p = (p x ,p y ). The BCS Hamiltonian 
(Al) acquires the following form in the real space 

& _ ( £(-*V - eA) -*(V • * + * • V)/2 \ 

l v -i(V-** + **- V)/2 _f(_iV + eA) J' 

(A2) 

It is now easy to check that the off-diagonal terms in 
Eq. (A2) are gauge-invariant under the simultaneous 
phase transformation of the fermion operators ip(r) — > 
i>(r) e %v ^ r \ ip^(r) — > tp^(r)e~ l<p ^ and the superconduct- 
ing order parameter M>(r) — > &(r) e 2tv ( r \ Notice that 
the gradient operators in Eq. (A2) act to the right on 
both ip(r) and \l/(r). 

Thus, it is incorrect to make the Peierls-Onsager sub- 
stitution p — > p — eA in the pairing potential A(p), as 
it was proposed in Refs. [40, 41]. Even if such a sub- 
stitution were made, it would have generated the term 
A ■ J d 3 pip(—p)vp(p) and its Hermitian conjugate. 
These terms vanish, because they reverse sign upon com- 
mutation of the fermion operators and changing sign of 
the variable of integration p. This conclusion also holds 
when the spin indices of the fermion operators ip are re- 
stored, because the spin-triplet order parameter * is a 
symmetric spin tensor, so there is no sign change upon 
exchange of the spin indices. Curiously, when the conven- 
tional vertex of interaction with the electromagnetic field 
is dressed by the impurity line as shown in Fig. 2b, the 
resulting vertices (14) and A y (15) have the structure 
similar to the terms discussed above. However, these ver- 
tices do not vanish upon fermion commutation, because 
they change sign upon exchange of the fermion frequen- 
cies wi and wi + ui n . 

The BCS Hamiltonian in the form (A2), with the pair- 
ing potential written as the anticommutator, was intro- 
duced and discussed in Ref. [7] and, more recently, in 
Rcf. [23]. This form was used in most papers on the 
subject, e.g., in Ref. [31], and was recently recognized 
in Ref. [44]. For a chiral d + id superconductor, a mi- 
croscopic derivation from a lattice model was given in 
Ref. [42]. 

The current operator is obtained by expanding the 
Hamiltonian (A2) to the first order in A, which gives 
the same result as in Eq. (6). Notice that the current op- 
erator for a superconducting systems is not given by the 
derivative dH /dp of the Hamiltonian (Al) with respect 
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to the momentum p, because p and A do not appear in 
the combination p — eA. For this reason, in contrast to 
non-superconducting systems, the anomalous Hall con- 
ductivity of a chiral superconductor is not expressed in 
terms of the Berry curvature, as discussed in the next 
Appendix. 



APPENDIX B: CHIRAL SUPERCONDUCTORS 
VS. THE TRSB TOPOLOGICAL INSULATORS 
AND METALS 

Insulators with a non-trivial band topology (the topo- 
logical insulators) can be divided into two classes: the 
time-reversal-invariant and the time-reversal-symmctry- 
breaking (the TRSB topological insulators) [73-76] . Here 
we discuss similarities and differences between chiral 
superconductors and the TRSB topological insulators. 
Some of their features, such as the energy gap in the 
single-particle spectrum and the presence of chiral edge 
states, are common. However, the electromagnetic prop- 
erties of chiral superconductors are different from those 
of topological insulators, e.g., because superconductors 
are not insulators. 

Following the pioneering work by Haldane [77], many 
authors investigated periodic systems with a topologi- 
cally non-trivial band structure resulting from a distribu- 
tion of the effective Aharonov-Bohm fluxes inside a unit 
cell. In such systems, the time-reversal symmetry may 
be broken even when the total flux through the unit cell 
is zero. Typically, such systems have a non-zero Berry 
curvature and exhibit anomalous Hall effect. A model of 
this kind, involving formation of the d xy + id x 2_ y 2 den- 
sity wave, was proposed in Rcf. [78] to explain the Kerr 
effect in underdoped cuprates [2]. The anomalous Hall 
effect was studied for this model earlier in Ref. [79], and 
the anomalous Nernst effect was discussed in Ref. [80]. 

However, it is important to emphasize that the origin 
of the anomalous Hall effect in these models has noth- 
ing to do with superconductivity. Thus, these models 
are not relevant for Sr 2 Ru04, where the Kerr effect ap- 
pears at the superconducting T c [1]. If the Kerr effect 
in Sr2Ru04 were due to a topologically non-trivial band 
structure, it would have been visible above the supercon- 
ducting T c , as in the underdoped cuprates [2]. Thus, for 
the applications to Sr2Ru04, the question is whether a 
non-zero anomalous Hall conductivity and the Kerr effect 
may originate from the TRSB solely due to chiral super- 
conductivity and not due to a topologically non-trivial 
band structure. For this purpose, one may consider the 
simple one-band parabolic dispersion law in the (x, y) 
plane and ignore complications of the real band structure 
of Sr 2 RuC>4, which consists of three distinct sheets [5, 9]. 
It was shown in Ref. [32] that the anomalous Hall con- 
ductivity of a clean p x +ip y superconductor vanishes even 
when an arbitrary electron dispersion is considered. 

For non-superconducting systems, the current opera- 
tor can be expressed in terms of the derivatives dH/dp 



of the Hamiltonian H with respect to the electron quasi- 
momenta p. For a topologically non-trivial band struc- 
ture, the off-diagonal matrix elements of the current 
operator between different bands produce a non-zero 
Berry curvature, and the anomalous Hall conductivity 
can be expressed completely in terms of the Berry cur- 
vature [46, 47, 77-80]. 

The situation is different for chiral superconductors. 
As discussed in Appendix A, the current operator for 
a superconductor cannot be expressed in terms of the 
derivatives dH/dp of the BCS Hamiltonian and is not re- 
lated to the group velocity dE(p)/dp of the Bogolyubov 
quasiparticles, as discussed around Eq. (3) in Ref. [43]. 
Thus, the Hall conductivity of a superconductor cannot 
be expressed in terms of the Berry curvature; moreover, 
the Hall conductivity vanishes, as discussed in Sec. I A. 
In contrast to the TRSB topological insulators, the effec- 
tive action of a chiral superconductor contains only one 
part of the Chern-Simons term [7, 26-32], involving Aq 
and A x or A y , but docs not contain the other part, in- 
volving A x and A y . The reason is that the superconduct- 
ing order parameter A has the phase which appears 
in the effective action and modifies requirements imposed 
by the gauge symmetry. When the self-consistent dynam- 
ics of the phase $ is taken into account, the anomalous 
Hall conductivity vanishes for a clean chiral supercon- 
ductor [26-29, 32, 37]. 



APPENDIX C: CALCULATION OF THE 
FUNCTION ff(w„) 

In the following Appendixes, we denote the gap A 
by simply A in order to shorten lengthy mathematical 
expressions. 

Here we calculate the function H{u n ) in Eq. (26): 



TUl 



\ 2 



yA 2 +^ v / A2 + (^+^) 2 / 

TUl + U n Wl 



TA J2 



[(^+ag 2 + A 2 ]3/ 2 ( ro 2 + A 2)3/2 



(CI) 



(C2) 



3n7; + uj n 



3tn; + 2u ri 



In going from Eq. (CI) to (C2), the terms odd in rui 
and vj\ + uo n vanish after summation over I. After the 
variable shift wi — > wi — uj n in the first and third terms 
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in Eq. (C2), H(cu n ) can be written as 



f(™l) 



Zwi — 2u) r , 



sjml + A 2 [(uji - iv, 
3n7; + 2uj n 



^fw'i + A 2 [{m 



A 2 } 



A 2 



(C3) 



(C4) 



The Matsubara sum in Eq. (C3) can be replaced by in- 
tegration in the complex plane 

r£/M-5 (C5) 

l J 

with the contour of integration being a series of circles 
around the points on the horizontal axis in Fig. 11, where 
the function tan(z/2T) has pole singularities. Then, this 
contour is expanded to a circle of an infinite radius and 
the contours encircling the poles zj and the brunch cuts 
of the function f(z), as shown in Fig. 11 

Tj2f(^) = \ E tan Res/O^) 

.3=1,2,3,4 



I 



/0)tan 



C 1 , O2 



/ z \ dz 
V2T1 ~2ni' 



(C6) 



The contribution from the poles at = —uj„ ± iA and 
Z 3A =w B ± iA is 

\ E tan (^) Re8 f^ = ^A taxih (w) (C7 ) 

j=l,2,3,4 



2A \2T r 
u> n + 3iA u> n — 3iA 



\ y/u) n (u n + 2iA) yjuj n (uj n - 2iA)^ 
The contribution from the branch cuts C\ and C 2 is 
If . . / z \ dz I f°° dx 

2] Cl J {z)i&n ^^ = -.] A 7^=^ 



(C8) 



3ix — 2w n 3ix + 2uj n \ l , ( x \ 



tanh 



A 2 + (ix - UJ n ) 2 A 2 + (ix + UJ n ) 2 J ~"""\2TJ 

Combining the two contributions (C7) and (C8) to 
Eq. (C3), we arrive to Eq. (27) for H(w n ). 



APPENDIX D: CALCULATION OF THE 
DRESSED VERTEX f d (p) 

Here we calculate the dressed vertex Tj (p) by solving 
the Bethe-Salpeter equation shown in Fig. 7b and given 
analytically in Eq. (45). According to Eq. (52), to cal- 
culate anomalous Hall response, we need to obtain the 
dressed vertex function T y . By iterating Eq. (45), we 
find a geometric series for T y 



fj/(p) = 7j/(p) + M y (zvi,w n ) 



(Dl) 




FIG. 11: (Color online) Contours of integration chosen to 
evaluate the Matsubara sums in Eqs. (C3) and (E5). The 
initial contour (blue lines) encircles the series of points along 
the horizontal axis, where tan(z/2T) has pole singularities. 
Then, the contour is deformed to infinity, encircling the poles 
Zj and the brunch cuts Ci and Ci of the function under the 
integral (red lines). 



where the function M y is 



M y (wi,u! n ) = rijUp hGp' (m)i y (p)Gp> (m + w„)r 3 
p' 

_ AyP F V F 



2-kt 



-L(wi, w n )(— ioifi + a 2 f 2 ). (D2) 



The functions ai, a 2 , and L are defined in Eqs. (47), (48), 
and (49), respectively. The next-order term in Eq. (Dl) 
is 

n i u l E hGp'(^i)M y (zui,uj n )g p/ (n7i + w„)f 3 
p' 

_ AyP F V F 



2-KT 



-L(mi,u) n ) (-iaih + a 2 f 2 ) (&if + ib 2 r 3 ), 



where the functions 61 and b 2 are defined in Eqs. (50) 
and (51). We observe that the higher-order terms bring 
the powers of the factor frifo + ib 2 f 3 , so we need to sum 
the following geometrical series 



g(Mb + *b 2 T,f = ^-LJ*_ T + {l _^ 2+hl 

(D3) 



T 3 . 



fe=0 



Using this result in Eq. (Dl), we find the dressed vertex 

r„(p) 

r y(p) = 1v(p) + ^ V ^ FVF - L(mi,u n )(-iaxTi + a 2 f 2 ) 
Zttt 

1 — bi ib 2 \ ,„ , 

x ! — -To + , XQ ; , 9 T 3 . (D4) 



(l-6!) 2 + 6 2 u (l-h) 2 + b 
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FIG. 12: (Color online) Contour of integration chosen to eval- 
uate the Matsubara sum in Eqs. (E2) and (F4). The initial 
contour (blue lines) encircles the series of points along the 
horizontal axis, where tan(z/2T) has pole singularities. Then, 
the contour is deformed to infinity and goes along the brunch 
cuts Ci,2,3,4 of the function under the integral (red lines). 



The contour of integration used to evaluate the sum in 
Eq. (E5) is shown in Fig. 11. The poles yield 



K 2a {uj n ) = itanh f-^j 



(E6) 



y/1 - (1 - iCo n f y/1 - (1 + iCj n ) 2 



\ 2i + Q n 
and the branch cuts contribute 



K 2 b K) 



-2 f — \/x 2 - ltanh 
Ji 2ir W 



2i — Cb n 



Ax 
2T 

1 



(E7) 



(2ix — w n )(l — [x + iidn} 2 ) (2ix+uj n )(l — [x — iuj n ] 2 ). 



Eqs. (E4), (E6), and (E7) give the expression for K(oj n ) 
shown in Eq. (61). 



APPENDIX E: CALCULATION OF THE 
FUNCTION K(u n ) 

The expression for the function K(uj n ) defined in 
Eq. (60) can be written as a sum of two terms 



K(u n ) 
Ki (w„) 



Kx{u) n ) + K 2 (uJ n ), 
TA 



(El) 



E 



2zui+uj n 



TA 

2tui + oj n 



v /(ro i +w„) 2 + A 2 7^ 



-A 2 
(E2) 

(E3) 



\J (ZUI + UJ n ) 2 + A 2 



[VJI + UJn) 2 + A 2 



One can notice that the above expressions for Ki(u> n ) 
and K 2 (oj n ) do not have a singularity at wi = — u> n /2. 
Similarly to Eq. (C3), the Matsubara sum in Eq. (E2) can 
be replaced by integration in the complex plane around 
the series of points along the horizontal axis in Fig. 12, 
where the function tan(z/2T) has pole singularities. The 
contour of integration is then deformed to go along the 
branch cuts 61,2,3,4 of the function under the integral, as 
shown in Fig. 12. The result is 

dx tanh(zA/2T) 



Ki(u n ) = -12cj n 



2tt ^ 



1 (w 2 + 4a; 2 ) 



(E4) 



where uj n = uj n / A. 

To evaluate the sum in Eq. (E3), we first shift the 
frequency vji — * zui + uj n in the first term to obtain: 



K 2 {u n )=TAY J 



tJvj 2 + A 2 



(2mi -u n ){[m -uj n ] 2 + A 2 ) 



y/wf + A 2 



(2wi + LU n ){[wi + LU n } 2 + A 2 ) 



(E5) 



APPENDIX F: CALCULATION OF THE 
FUNCTION J(oj n ) 

The Matsubara sum (71) for the function J(ui n ) can 
be separated into the terms with poles and branch cuts 



jK) = taW * 



1 



zu 2 + A 2 (mi + uj n ) 2 + A 2 



(Fl) 



^{m l +Lu n ) 2 + A 2 ^zu 2 + A 2 ) ' 
After shifting wi in the second term of Eq. (Fl), we get 
2TA 2TA \ 



( 2' 



+ A 2 ^/(vu l +oj n ) 2 + A 2 y/w 2 + A 2 J ■ 

(F2) 



The first term in Eq. (F2) has only simple poles 
^ — ft—: o = tanh ( | . 



wf + A 2 



2T 



(F3) 



The contribution from the second term in Eq. (F2) can 
be calculated by using the contour shown in Fig. 12 



E 



2TA 



— / \ 

2 /^tanh(^) 



^ ^{vj l +uj rL ) 2 + A 2 ^vo 2 + A 2 it J V* 2 - A 2 

A 



A 



A 



y/(ix + u) n ) 2 + A 2 sjiix - uj n ) 2 + A 2 



(F4) 



By combining the two contributions (F3) and (F4), we 
obtain the result shown in Eq. (72). 
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